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The complete 2-loop quantum electrodynamic corrections to the muon lifetime are calculated 
in the Fermi theory. The exact result for the effects of virtual and real photons, virtual electrons, 
muons as well as e''"e“ pair creation is 
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where To is the tree-level width. The theoretical error in the value of the Fermi coupling con¬ 
stant, Gf, is now rendered negligible compared to the experimental uncertainty coming from the 
measurement of the muon lifetime. The overall error in Gf is then roughly halved giving 


Gf = (1.16637 ±0.00001) x 10“® GeV^ 
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I. INTRODUCTION 

The Fermi coupling constant, Gp, plays a key role in precision tests of the Standard Model of electroweak interac¬ 
tions. Along with the electromagnetic coupling constant, a, and the Z boson mass, Mz, it is one of the best measured 
quantities of electroweak physics and as such is used as input in all higher-order calculations. Gf is one of the few 
quantities that is sensitive to physics at very high energy scales and is intimately related to the p parameter . It was 
the value of Gf that provided some of the strongest constraints on the mass of the top quark before it was directly 
observed. 

Gf is extracted from measurements of the muon lifetime, = F”^, which is a purely leptonic process and therefore 
very clean both experimentally and theoretically. Its quoted error is SGp/Gp = 1.7 x 10“^ of which 0.9 x 10“® is 
experimental and 1.5x10“® is theoretical; the latter being an estimate of the size of the 2-loop corrections. Experiments 
are under consideration at Brookhaven National Laboratory, the Paul Scherrer Institute, and the Rutherford-Appleton 
Laboratory that could lead to a reduction in the experimental error on the of a factor of 10 or more. 

The radiative corrections to muon decay in the full Standard Model naturally factorize into two pieces |^, one of 
which, to a very high degree accuracy, is just the quantum electrodynamic (QED) radiative corrections in the Fermi 
theory. The other piece is left free of infrared singular contributions. It contains purely weak corrections that can be 
absorbed into Gp which then possesses an enriched sample of weak sector physics. Such a separation between QED 
and weak corrections is not generally possible for charged current processes. 

The 1-loop QED contributions to the muon lifetime were first calculated over 40 years ago by Kinoshita and Sirlin 
and by Berman [Q. It is known Q that the Fermi theory in the presence of QED is finite to leading order in Gf 
and to all orders in the electromagnetic coupling constant, a. This remarkable fact means that Gp can be defined in 
a physically unambiguous manner at least up to the point where finite W propagator effects begin to appear. 

In this article the 2-loop QED radiative corrections to muon lifetime are calculated in the Fermi theory. The result 
is used to extract an improved value for Gf in which the error is entirely due to the experimental uncertainty. 


II. THE FERMI COUPLING CONSTANT 

The Fermi theory Lagrangian, relevant for the calculation of the muon lifetime is 
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Here Cw is the Fermi contact interaction 

Cw = -2V2 Gf ■ [^e7A7LV’°J (2) 

in which '0^, 'ijjf., and are the wave functions for the muon, the electron and their associated neutrinos 
respectively. The Euclidean metric in which time-like momenta squared are negative is used. /Iqcd i® bare 
Quantum Chromodynamic (QCD) Lagrangian responsible for the strong interactions and >Cqed i® usual bare 
Lagrangian of QED, 

-CgED = - X! ^ {dpAl - daA^pf -f ^ Qfp^°fjpp)°fA°. (3) 

/ / 

The sum is over all fermion species, /, with mass, rrif, and electric charge, Qf. Ap is the photon field and 7 l = 5 ( 1 - 1 - 75 ) 
denotes the usual Dirac left-hand projection operator. The superscript ° indicates bare, as opposed to renormalized, 
quantities. For the present purposes Gp goes unrenormalized. Throughout this article dimensional regularization 
is used for the ultraviolet (UV) divergences. The appearance of infrared (IR) divergences is largely avoided by the 
methods employed here. 

The formula obtained for Tp by means of the Cp is finite to leading order in Gp and all orders in the renormalized 
electromagnetic coupling constant, ar = e^/(47r) Q. This follows from the fact that under a Fierz rearrangement that 
interchanges the wavefunctions tpe and tpi,^ in Cw the currents remain purely left-handed vector currents. This is in 
sharp contrast to the case of neutron decay in which scalar and pseudoscalar terms are generated and for which the 
following arguments break down. The radiative corrections in that case are not hnite. Considering the vector part, 
ipelp'ipp, of this effective /i-e current, one sees that after fermion mass renormalization is performed the remaining 
divergences are independent of the masses and thus cancel, as for the case of pure QED. The QED corrections to the 
axial vector part may be shown to be finite by noting that the transformations tpe —>■ qsV'e and nie —>■ —me leave 
'Cqed and £qcd invariant but exchanges tpejxtpp <-^ tpejxj^tpp. Thus the radiative corrections to the axial-vector 
part of the current are equal to those of the vector part in the limit of me = 0. In practice then the only radiative 
corrections to the vector pieces in Cw need to be calculated which avoids entirely the problems associated with 75 in 
dimensional regularization. 

To lowest order in Ge the expression for the muon lifetime calculated from Cp takes the general form 

-=r^ = ro(l + Ag). (4) 
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and Ag encapsulates the higher order QED and QCD corrections generated by £e and can be expressed as a power 
series expansion in the renormalized electromagnetic coupling constant, Or, 

00 

A<7 = ^Ag« (5) 

2 = 0 


in which the index i gives the power of olt that appears in Agb). 

Assuming that the electron neutrino and muon neutrino are massless it can be shown that 
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that comes from phase space integrations. 

The 0{a) corrections in Ag, first obtained by Kinoshita and Sirlin H and by Berman 0, are 
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where ^ is the Riemann zeta function and (^(2) = An exact expression for the full electron mass dependence in 

has been given by Nir 0. 

Recently the hadronic contributions to were computed using dispersion relations along with contributions 

from muon and tau loops 1^. Their effect was shown to be small relative to the present experimental error. They 
become relevant for the next generation of muon lifetime experiments but the hadronic uncertainty is still well under 
control. 

The Kinoshita-Lee-Nauenberg ^ theorem guarantees that Aq is free from singularities as rrie —> 0, other than those 
that can are absorbed into ar- It may be shown that all large logarithms of the form Q;Mn*“^(m^/TOg) for all 
i > 0 and those of ln{m^/rnl) can be accounted for, in a manner consistent both with the calculation of Ref. 
and the perturbative results presented here, by setting 
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where a is the experimentally-measured quantity, a = 1/137.0359895(61) |lT| . The contribution to the muon lifetime 
from the 0{a^) logarithmic term coincides with the result obtained in Ref. |^. The logarithms of 0[a^) were first 
obtained by dost and Luttinger jl^]. When evaluated Eq/I) yields aeinifi) = 1/135.90 = 0.0073582. In the MS 
renormalization scheme with ’t Hooft mass, /r = Eq.(p) correctly includes non-logarithmic terms up to 0{a^) 
but those of 0{a^) have been dropped. 


III. 2-LOOP CORRECTIONS 


A. Photonic Corrections 

The calculation of the 2-loop QED corrections to the muon lifetime involves the sum of the cross-sections ^ 
v^e~Ve, fJ.~ —> r'^e“7'e7) ^ v^e~Vell and iJ.~ —> Vi^e~Vf.e'^e~ with up to two virtual photons. Individual diagrams 
are IR divergent and in some cases require integration over a 5-body phase space. The problem of cancelling these IR 
singularities can be avoided entirely if the QED corrections are obtained as the imaginary part of 4-loop propagator 
type Feynman diagrams by means of the optical theorem. Some of these diagrams are shown in Fig.l. The heavy lines 
represent muons which are the only particles taken to have non-zero mass. The 4-fermion vertex used is the vector 
part of the usual one from the Fermi theory. Inspection of the diagrams shows that the cuts generating imaginary 
parts produce all of the Feynman diagrams contributing to muon decay. Extra diagrams do appear in which the cut 
goes through a muon line but such diagrams vanish kinematically because the external muon is on its mass shell. 




(a) 


(b) 





FIG. 1. Examples of diagrams whose cuts give contributions to p —> Ve, M ^ or p ^ i^e77. 

The imaginary parts of the necessary 4-loop propagator type diagrams were calculated as follows. Recursion relations 
Ip^ obtained by integration-by-parts were first applied to reduce all dimensionally regularized integrals to a small set 
of relatively simple integrals. These primitive integrals were chosen to be free from specific IR divergences that occur 
on-shell. The well-behaved primitive integrals were then calculated by taking the external muon momentum, q, off 
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mass shell to obtain expressions as power series in a; = —q^/rn^ and logarithms of x using well-established large mass 
expansion techniques along the lines of Ref. . This series serves as a convenient representation as its coefficients 
involve simpler integrals. Many terms in the large mass expansion can be discarded since they do not contribute to 
the imaginary part. What remains of the coefhcients in the expansion can be evaluated in closed form in terms of 
polygamma functions and certain classes of multiple nested sums fl^] . Then, the on-shell limit, a; = 1, is taken and 
the infinite sum over the coefficients of x^ is performed. In this process the exact expressions collapse into known 
constants such as the Riemann zeta function of integer arguments, and polylogarithms, Life(l/2). Details of the 
procedures followed will be described elsewhere . 

Fermion mass renormalization is performed in the on-shell scheme]^ that generates derivatives of fermion self-energies 
for the external leg corrections. All diagrams were calculated in a general covariant gauge for the photon held and 
exact cancellation in the hnal result of the dependence on the gauge parameter was demonstrated. 

The result for just the photonic diagrams is 




= rn 
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B. Electron Loop Corrections 




FIG. 
or fj,~ - 


(c) (d) 

2. Diagrams containing an electron loop whose cuts give contributions to muon decay, 
-+ u^e~ e~. 


Ufj,e Ve, /r 


v^l,e Ue-y 


The contribution of electron loops to the muon lifetime differs from those of other fermions in that they must be 
combined with diagrams with an additional e+e“ pair in the hnal state in order to produce an IR hnite result, however 
the procedure described above may be applied here as well. The electron loop diagrams are shown in Fig.|^. To Fig.^jd 
must be added a diagram containing a muon mass counterterm, 6m^, on the external leg. Furthermore, diagrams, in 
which the electron loop is replaced by the photon self-energy counterterm, must be included to produce a UV hnite 
result. This counterterm contribution is proportional to and depends on the particular renormalization scheme 

that has been chosen. The overall result in the MS renormalization scheme with ’t Hooft mass p, = consistent 
with Eq.(||) is 
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^That is to say that the renormalized mass of a stable fermion is set equal to its physical or pole mass 
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whi ch is about two orders of magnitude greater than that of either muon loops or hadrons. The value obtained in 
Eg.( 10b ) is consistent with a numerical study presented in Ref. 0 in the context of semi-leptonic decays of heavy 
quarks. 

The same methods used to calculate the contribution from electron loops can be applied to muon loops. Agreement 
was found with the result of Ref. [HI. 


IV. CONCLUSIONS 


The photonic corrections of section III A can be combined with those of the electron loops and e“'"e pair production 
of section III B, and adding the exact result for muon loops of Ref. gives 


AT 


( 2 ) 

QED 


/a^\2/156815 
V 5184 

= ro (^)%.743 


1036 

27 


C(2)-^C(3) + yC(4) + 53C(2)ln(2)) 


(lla) 

(llb) 


with ar = ae(TO/j) = l/135.90. The resulting expression contains all corrections of C>(a^), O (a^ ln(mg/m^)) and 
O (aMn*“^(TOg/TO^)) for all i > 2. Adding the hadronic and tau loop contributions of Ref. one obtains 


Ar(2) = To ^ (6.700 ± 0.002). 


( 12 ) 


where the error is a conservative estimate of the hadronic uncertainty. Using the current best value for = (2.19703± 
0.00004) /rs 1^ yields 


Gf = (1.16637 ±0.00001) x 10“® GeV^ 


(13) 


That represents a reduction in the overall error on Gf of about a factor of 2 and a downward shift in the central 
value of twice the experimental uncertainty. Gp is now known to 9 ppm. The next generation of measurements of the 
muon lifetime are expected to reduce this by at least a further factor of 10. 
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